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Doppler  shift  measurements  of  the  Na  Di  absorption  line  have  revealed  solar  oscil¬ 
lations  in  a  new  regime  of  fre(iuency  and  wavenumber.  Oscillations  of  vertical  velocities 
in  the  temperature  minimum  and  low  chromosphere  of  the  Sun  are  observed  with  fre¬ 
quencies  ranging  up  to  9. -5  mHz.  The  fundamental  modes  appear  with  wavenumbers 
up  to  5.33  Mm~^  (equivalent  spherical  harmonic  degree,  3710).  We  find  no  evidence 
for  chromospheric  modes  of  3  minute  period. 
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Doppler  shifts  and  intensities  were  measured  for  the  core  of  the  Na  D;  absorption  line 
using  the  Orbiting  Solar  Laboratoiw  evaluation, -tunable-filter  system,  and  CCD  camera 
(Title  et  al.  1991)  at  the  Swedish  Solar  Observatory.  Filtered  images  (filtergrams)  were  taken 
at  three  wavelength  offsets  (—85.0, +80  m.\)  from  the  nominal  line  center.  The  filter  width 
was  ~  62  m.4  (FWHM).  Hence,  the  measurement  was  sensitive  to  line-of-sight  velocities 
500-900  km  above  the  photosphere  (Schleicher  1976). 


Active  region  AR  6181  was  observed  for  2.4  hr  on  1990  .August  2  using  the  Swedish  50 
cm  vacuum  solar  telescope  on  La  Palma,  Spain.  The  224  ’  x  224  ’  field  of  view  was  resolved 
into  512  X  512  pixels.  .A  set  of  filtergrams  was  taken  <'ver\  17.7  s  for  a  total  of  480  pairs  of 
Doppler  measurements  ( Dopplergrains)  and  intensity  measurements.  The  angular  size  of  the 
field  of  view  was  determined  to  within  1%  by  comparison  of  magnetograms  from  La  Palma 
and  from  the  National  Solar  Observatory  at  Kitt  Peak. 


Reduction  of  filtergrams  included  calibration  for  the  CCD  detector  response  and  removal 
of  image  rotation  and  drift  caused  by  the  telescope  tracking.  The  effect  of  atmospheric 
“seeing”  was  minimized  by  “de-stretching”  (November  1986)  to  remove  image  distortion,  and 
spatial  filtering  to  normalize  the  modulation  transfer  function  (MTF)  of  the  atmosphere.  The 
“MTF  normalization"  forces  the  spatial  spectrum  of  each  filtergram  of  a  given  wavelength 
to  look  like  the  average  spectrum  for  filtergrams  of  that  wavelength.  Doppler  shifts  and  line- 
center  intensities  were  then  calculated  at  each  pixel  by  solving  for  the  vertex  of  a  parabola 
through  the  three  filtergram  intensities. 

Three-dimensional  Fourier  transforms  were  made  of  the  sets  of  Dopplergrams  and  inten¬ 
sities  after  apodization  in  the  space  and  time  directions.  Both  the  full  field  of  view  and 
a  79”  X  79”  subfield  of  "ciuiet"  Sun  were  transformed.  The  temporal  and  spatial  Nyquist 
frequencies  were  28.3  mHz  and  9.9  Mm~*,  respectively.  The  distribution  of  power  in  Fourier 
{k,u)  space  as  expected  for  acoustic  waves  (p-modes)  and  surface  waves  (/-modes)  defines 
“trumpet”-shaped  surfaces  (Hill  1988)  of  approximate  form  u  oc  ||A"||2. 

Because  the  telescope  was  tracking  the  active  region,  there  is  no  net  motion  of  the  medium 
across  the  field  of  view.  Hence,  the  frequency  of  a  “trumpet"  at  wavenumber  k  is  the  same 
as  the  frequency  at  ~k.  The  active  region,  however,  was  17°  from  disk-center,  causing 
waves  propagating  toward  or  away  from  disk  center  to  appear  at  higher  wavenumber  than 
waves  propagating  in  the  perpendicular  direction.  This  foreshortening  distorts  the  cylindrical 
symmetry  of  the  “trumpets"  to  give  them  an  elliptical  rather  than  circular  cross  section. 
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The  power  spectra  were  integrated  in  tlie  azimuthal  direction  with  allowance  made  for 
the  foreshortening.  The  two-dimensional  k'-i-/  spectrum  for  the  full-field  Dopplergrams  is 
shown  in  Figure  1.  The  p-mode  and  /-mode  ridges  are  clearly  visible  at  5.5  mHz,  but  there 
is  not  a  horizontal  band  of  jjowei  nc"  are  there  "avoided  crossings”  as  would  be  expected  for 
a  3  minute  chromospheric  resonance  (Ulrich  Sz  Rhodes  1977). 

In  addition  to  producing  the  power  spectrum,  the  three-dimensional  Fourier  transform 
wa5  used  to  study  the  phase  of  the  oscillations.  The  pha.se  difference  between  velocity  and 
intensity  oscillations  was  measured  to  be  ~  0°  for  ail  visible  ridges. 

Slices  of  the  power  spectra  were  fitted  using  a  nonlinear,  weighted,  least-squares  fitting 
program.  The  weights  were  determined  by  the  statistical  error  on  the  measured  power, 
which  in  turn  was  derived  from  an  assumed  error  of  50  m/s  in  the  velocity  measurement. 
The  assumed  power  spectrum  shape  consisted  of  a  power-law  background,  Pkiiz)  = 
together  with  Gaussian  ridges. 

The  I  kground  paiameters.  and  were  fitted  independently  for  each  wavenumber 
slice  of  h  spectrum.  The  fit  included  a  range  of  frequencies  on  either  side  of  the  ridges: 
(0.12  <  fc'/mHz  <  1.3)  and  (12.0  <  ujmWz  <  15.0).  The  background  velocity  power, 
integrated  over  frequency,  is  plotted  in  Figure  2.  The  spectrum  falls  off  much  more  slowly 
with  wavenumber  than  the  spectrum  of  convection  in  the  photosphere  reported  by  Chou  et 
al.  (1991).  We  find  a  root  mean  square  velocity  of  about  820  m/s  for  convection  in  the  low 
chromosphere. 

The  ridges  were  fitted  simultaneously  after  subtracting  the  background.  Frequencies  of 
ridges  from  the  velocity  spectra  are  plotted  in  Figure  3  with  their  respective  statistical  errors. 
All  points  plotted  correspond  to  fits  in  which  the  offset,  width,  and  height  of  the  Gaussian 
converged  to  physically  reasonable  values,  and  the  integrated  power  under  the  Gaussian  was 
at  least  4  standard  deviations  above  the  background. 

The  mode  power  is  plotted  as  a  function  of  frecpiency  in  Figure  4.  Mode  power  is 
defined  as  the  integral  under  the  fitted  Gaussian.  The  p-mode  power  has  the  same  frequency 
dependence  for  all  p-ridges.  There  is  very  little  explicit  wavenumber  dependence.  Power 
increases  as  to  a  maximum  at  3.4  mHz,  at  which  point  it  turns  over  and  falls  off  again 
as  The  /-modes  have  less  power  at  a  given  frequency  than  the  p-modes.  As  surface 

gravity  waves  with  very  little  compression  they  are  physically  different  from  the  p-modes 
and  are  not  expected  to  show  the  same  power  spectrum. 

An  MTF  correction  was  made  to  the  power  plotted  in  Figures  2  and  4  by  combining  the 
calculated  MTF  of  the  telescope  and  CCD  with  an  estimated  atmospheric  MTF.  The  latter 
was  derived  from  the  difference  betwet^n  power  spectra  of  the  ".M  IT'  normalized”  filtergrams 
and  the  sharpest  filtergrams  in  the  set.  Since  even  the  best  images  will  still  include  some 
atmospheric  blurring,  this  correc  tion  will  provide  an  underestimate  for  power  at  the  highest 
wavenumbers. 


All  six  fitted  p-mode  ridges  c-xtead  up  to  at  least  7.0  mHz.  The  ps  and  pe  ridges  extend 
up  to  a  frequency  of  9.5  mHz.  We  see  no  deviation  from  a  smooth  extrapolation  of  the 
ridges  at  lower  frequencies.  This  suggests  that  there  are  no  new  cavity  resonances  in  the 
high  frequency  range  revealed  in  this  work.  Libbrecht,  Woodard  &  Kaufman  (1990)  have 
made  measurements  that  o\’erlap  with  ours  at  frequencies  up  to  5.3  mHz.  The  two  data  sets 
agree  to  within  the  experimental  scatter  of  ±15  //Hz. 

The  /-ridge  is  \  isible  at  wa\enuml)ers  ii]>  to  5.33  Min~’.  There  is.  however,  a  systematic 
discrepancy  between  fre(|uencies  filled  from  the  velocity  and  intensity  spectra.  Frequencies 
measured  from  the  velocity  spectra  are  consistently  higher  than  those' from  the  intensity 
spectra  by  up  to  40  //Hz  at  wa\enumbers  below  2.5  Mm“^  This  difference  increases  to  ~ 
100  //Hz  at  4  Mm~'.  The  discrepancy  may  in  part  be  due  to  a  difference  between  background 
shapes  in  the  velocity  and  intensity  spectra. 

Libbrecht  tt  al.  (1990)  report  a  Gaussian  background,  centered  near  3  mHz,  in  their 
velocity  spectrum.  Inclusion  of  such  a  background  in  our  fit  would  bring  the /-ridge  frequen¬ 
cies  into  better  agreement  at  wavenumbers  below  1.1  Mm“*.  but  would  push  them  farther 
apart  at  higher  wavenumbers.  Fitting  such  a  background  function,  moreover,  would  require 
a  longer  time  series  with  correspondingly  better  signal-to-noise  ratio.  It  is  noteworthy  that 
the  /-mode  frequencies  measured  from  the  intensity  spectra  agree  with  those  reported  by 
Libbrecht  e/  al.  to  within  15  //Hz. 

The  discrepancy  could  also  be  caused  by  a  difference  in  mode  frequencies  in  the  magnetic 
and  nonmagnetic  regions  within  the  field  of  view.  The  weighting  of  magnetic  regions  is 
different  for  intensity  and  velocity  s|)ectra. 

Regardless  of  which  set  of  /-mode  frequencies  we  use.  the  measured  frequencies  lie  below 
the  theoretical  frequencies  given  by  =  gk  at  wavenumbers  greater  than  1  Mm“'.  The 
difference  is  100-200  //Hz  at  4  Mm“'.  The  source  of  the  frequency  shift  is  unknown.  Camp¬ 
bell  Roberts  (1989)  and  Evans  Robeits  (1990,  1991)  predict  that  a  horizontal  magnetic 
field  in  the  chromosphere  w’ould  increa.se  the  frequencies  of  the  /-modes  by  increasing  the 
restoring  force  on  the  surface  waves.  The  sign  of  the  frequency  shift  may  reverse  at  high 
wavenumber,  however,  in  a  mor<*  realistic  model  of  the  chromosphere  with  an  appropriate 
magnetic  field  structure  (Roberts  1991). 

Duvall  et  al.  ( 1991 )  report  frequencies  up  to  6.7  mHz  for  //-modes  in  a  lower  wavenumber 
regime  than  we  can  explore  with  our  data.  Their  measurements  extend  above  the  acoustic 
cutoff  frequency  of  ~  5.3  mHz.  At  high  wavenumber  the  acoustic  cutoff  frequency  is  ap¬ 
proximately  proportional  to  waienumber  (Gough  1990).  Hence,  our  measured  frequencies 
all  lie  below  the  acoustic  cutoff.  This  is  borne  out  by  the  measured  0°  phase  difference 
between  velocity  and  intensity  oscillations.  This  phase  difference  is  indicative  of  standing 
waves  (Deubner  <1:  Fleck  1989). 

Two  explanations  have  been  put  forward  for  the  origin  of  the  high  frequency  ridges  ob¬ 
served  by  Duvall  et  al.  (1991).  Kumar  et  al.  (1990)  propose  that  waves  are  emitted  isotrop- 


ically  from  a  source  near  the  tO])  of  the  convection  zone.  Waves  traveling  directly  toward 
an  overlying  observed  surface  interfere  with  waves  reflected  from  below  to  produce  “mock 
modes.”  Balmforth  k  Gough  (1990),  on  the  other  hand,  propose  a  chromospheric  cavity 
which  couples  to  (he  familiar  subphotospheric  cavity  to  allow  high  frequency  resonances. 

The  clear  absence  of  chromosiiheric  resonances  in  our  data  casts  doubt  on  the  existence  of 
the  chromospheric  cavity  required  by  the  model  of  Balmforth  k  Gough.  Either  the  reflection 
coefficient  at  the  ba.se  of  the  corona  is  too  low  or  the  surface  is  too  irregular  to  form  a  good 
cavity. 
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Figure  1:  Scaled  velocity-power  spectrum.  \  portion  of  the  collapsed  power  spectrum  for 
the  full  field  of  view,  with  variable  gray-scaling  in  the  A-  and  directions.  Background  has 
been  subtracted.  The  noise  visible  at  the  top  of  the  spectrum  is  an  artifact  of  the  scaling. 


Figure  2:  Background  \elocity  power  integrated  over  frecpiency.  Background  was  fitted  from 
the  “quiet”  subfield  spectrum.  S<)me  correction  has  been  made  for  the  MTF  of  the  telescope 
and  seeing  (see  text). 


Figure  3:  Measured  mode  frequencies  (from  velocit}'  spectra):  //  vs.  logjp  A*.  Frequencies  of 
the  / and  pi~pe  ridges  are  plotted  with  1  a  error  bars.  For  k  <  3.96  Mm“’,  frequencies  were 
fitted  from  the  full-field  sjK'ctrum.  At  highcM-  k  the  "(|ui('t‘‘  subfield  spectrum  was  used,  as 
it  had  less  noise. 


Figure  4:  Velocity  power  \  s.  mode  frecnieiicy.  Some  correction  has  been  made  for  the  MTF 
of  the  telescope  and  seeing  (see  text). 
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